Interpolation engine for analysis of time-varying load data signals

ABSTRACT

A method for analyzing fatigue life of an elastomeric component includes a step of conducting a finite element analysis to obtain a base state. A plurality of case vectors are then selected to represent a space of possible loading states that occur within a time-varying load data signal based on measurement of the elastomeric component or vehicle dynamics. For at least a portion of the case vectors, a finite element analysis is conducted at a plurality of discrete gridpoints along the case vectors starting at the base state and tracking the case vector. Using an interpolation engine, desired local solution variables for a current state may be interpolated from the finite element analysis at the plurality of discrete gridpoints. A damage calculation may then be calculated based on the desired local solution variables for the current state.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. Utility application Ser. No.13/760,169, filed on Feb. 6, 2013, and issued as U.S. Pat. No. 9,645,041on May 9, 2017, which claims the benefit of U.S. Provisional ApplicationNo. 61/595,329, filed on Feb. 6, 2012. The entire disclosures of theabove applications are hereby incorporated herein by reference.

FIELD OF THE INVENTION

The present disclosure relates to a method and system for analyzingfatigue life of elastomeric or rubber components.

BACKGROUND OF THE INVENTION

Solutions for fatigue analysis from finite element analysis (FEA) ofmetallic components have been available for many years. An importantpart of the analysis for linear structures such as metallic componentsis a procedure called “scale and combine”, which allows one to convertfrom raw road load data to stresses and strains. In this procedure, aseries of unit load cases is modeled in FEA and can then be used toreconstruct stress or strain histories for a multiaxial input signal.

Nonlimiting examples of fatigue analysis solutions for metalliccomponents are described in each of: Conle, F. A., and C-C. Chu.“Fatigue analysis and the local stress-strain approach in complexvehicular structures.” International journal of fatigue 19.93 (1997):317-323; Braschel, Reinhold, Manfred Miksch, and Rolf Schiffer. “Methodof monitoring fatigue of structural component parts, for example, innuclear power plants.” U.S. Pat. No. 4,764,882. 16 Aug. 1988; Yim, HongJae, and Sang Beom Lee. “An integrated CAE system for dynamic stress andfatigue life prediction of mechanical systems.” Journal of MechanicalScience and Technology 10.2 (1996): 158-168; and Conle, F. A., and C. W.Mousseau. “Using vehicle dynamics simulations and finite-element resultsto generate fatigue life contours for chassis components.” Internationaljournal of fatigue 13.3 (1991): 195-205.

For many, FEA has become an essential part of maturing and qualifyingdesign concepts, providing a cost-effective and proven basis forjustifying investment in physical prototypes and testing. However,conventional fatigue analysis solutions do not work well for elastomericcomponents because of their macromolecular structure. In particular, thescale and combine method is not suitable for rubber parts, because ofmaterial and kinematic nonlinearities in rubber.

Rubber or elastomeric components exhibit unique behavior and requirespecialized analysis methods. Developing a durable elastomeric componentoften involves expensive, time-consuming, trial-and-error iterations.There has been a long-felt, but unsolved, need in industries such asautomotive, defense, transportation, heavy equipment, offshore, medicaldevices and consumer products, for a solution to put developers incontrol of durability issues early in the development cycle, when thegreatest opportunities to influence performance exists.

With regard to rubber components such as bushings, tire treads, seals,etc. used in an automotive setting, it is known that road load signalsare too lengthy to use for full FEA of the rubber components. However, afull strain history remains desirable for damage calculations by FEA.

There is a continuing need for a method and system for efficientlyobtaining strain and stress histories at potential failure locations ina rubber component, based on a given time-varying load data signal suchas a road load input signal and FEA.

SUMMARY OF THE INVENTION

In concordance with the instant disclosure, a method and system forefficiently obtaining strain and stress histories at potential failurelocations in a rubber component, based on a given a time-varying loaddata signal such as a road load input signal and FEA, is surprisinglydiscovered.

In one embodiment, a method for analyzing fatigue life of an elastomericcomponent includes analyzing a time-varying load data signal obtained,e.g., from measurement of the elastomeric component or from a multibodydynamics analysis of the elastomeric component. A finite elementanalysis of the elastomeric component is conducted to obtain a basestate. A plurality of case vectors are then selected to represent aspace of possible loading states that occur within the time-varying loaddata signal. For at least a portion of the case vectors, a finiteelement analysis is conducted at a plurality of discrete gridpointsalong the case vectors, the gridpoints selected along the case vectorsstarting at the base state and tracking the case vector. Using aninterpolation engine, desired local solution variables for a currentstate may be interpolated from the finite element analysis at theplurality of discrete gridpoints. A damage calculation may then beperformed based on the desired local solution variables for the currentstate.

In another embodiment, a method for analyzing fatigue life by analysisof a time-varying load data signal obtained from measurement of anelastomeric component includes a step of identifying independentvariables of the time-varying load data signal. A finite elementanalysis is then conducted to obtain a base state of the elastomericcomponent. A plurality of case vectors are then selected to represent aspace of possible loading states that occur within the time-varying loaddata signal. For each of the case vectors, a finite element analysis isalso conducted at a plurality of gridpoints along the case vectors, thegridpoints selected starting at the base state and tracking the casevectors. Using an interpolation engine, desired local solution variablesfor a desired current state are then obtained. The local solutionvariables provide at least one of an interpolated strain history and aninterpolated stress history. A damage calculation is then performedbased on the one of the interpolated strain history and the interpolatedstress history, in order to determine a potential failure location inthe elastomeric component.

In a further embodiment, a system for analysis of a time-varying loaddata signal includes an interpolation engine. The interpolation enginehas at least one processor and at least one memory. The at least onememory includes a computer readable medium having a set ofcomputer-readable instructions embodied thereon that, when executed bythe at least one processor, cause the at least one processor to performa method according to the present disclosure. In particular embodiments,the processor performs a method of: interpolating a strain history and astress history of the elastomeric component at a current state from asimplex defined by neighboring case vectors radiating outwardly from abase state of the elastomeric component. The case vectors represent aspace of possible loading states that occur within a time-varying datasignal obtained from measurement of loads on the elastomeric component.The case vectors also have a plurality of discrete gridpoints disposedthereon. The interpolated strain history and the interpolated stresshistory may then be used in performing a damage calculation to determinefatigue life and a potential failure location in the elastomericcomponent.

DRAWINGS

The above, as well as other advantages of the present disclosure, willbecome readily apparent to those skilled in the art from the followingdetailed description, particularly when considered in the light of thedrawings described herein.

FIG. 1 is a graphical illustration of an exemplary dual channeltime-varying load data signal for a rubber component undergoing a loadcycle, the time-varying load data signal for use with the method andsystem of the present disclosure;

FIG. 2 is a graphical illustration of the dual channel time-varying loaddata signal shown in FIG. 1, with a first channel disposed along anX-axis and a second channel disposed along a Y-axis;

FIG. 3 is a graphical illustration of the method of the presentdisclosure, showing a selected base state, a series of case vectors, anda plurality of gridpoints along the case vectors which are used for aninterpolation of a current state of the rubber component;

FIG. 4 is a flow chart illustration showing a method for analyzingfatigue life of a rubber component, according to one embodiment of thepresent disclosure;

FIG. 5 is a flow chart illustration showing a further method forinterpolating of strain and stress histories as part of the method shownin FIG. 4; and

FIG. 6 is a schematic illustration of a system for analyzing fatiguelife of a rubber component, according to one embodiment of the presentdisclosure.

DETAILED DESCRIPTION OF THE INVENTION

The following detailed description and appended drawings describe andillustrate various embodiments of the invention. The description anddrawings serve to enable one skilled in the art to make and use theinvention, and are not intended to limit the scope of the invention inany manner. In respect of the methods disclosed, the order of the stepspresented is exemplary in nature, and thus, is not necessary or criticalunless otherwise disclosed.

The term “road load”, as is used herein, applies to any long load ordisplacement signal containing random content or containing varyingfrequency, amplitude, and phase. For example, the road load mightoriginate from measurements made on the road, in the air, or in any typeof service where there are time-varying loads on a component.

A method for analyzing road load data to select a subset forcomputer-aided engineering (CAE) analysis is described in U.S. Pat.Appl. Publication No. 2004/0254772 to Su, the entire disclosure of whichis hereby incorporated herein by reference. A method for computing roadload history from a vehicle dynamics model, using a particular approachfor modeling tire behavior, is also described in U.S. Pat. No. 7,363,805to Jayakumar et al., the entire disclosure of which is herebyincorporated herein by reference. Other suitable methods for acquiringand processing road load data or other time-varying load data may alsobe used within the scope of the present disclosure.

An example is made herein of an elastomeric or rubber component in theform of a simple rubber bushing in an automotive context, undergoing thetime varying load in the form of the road load through operation of avehicle having the rubber bushing. However, it should be understood thatthe method and system of the invention may be used to predict potentialfailure locations of any rubber component to which a time-varying loadis applied in service, for example, a tire component, such as a rubbertread, an engine mount, a rubber seal, a rubber track, etc. Othersuitable types of rubber components may also be analyzed for potentialfailure modes and locations using the method and system of the presentdisclosure.

With reference to FIGS. 1-6, a unique method and system for efficientlyobtaining strain and stress histories at potential failure locations ina rubber component, such as a rubber bushing, based on a given atime-varying load data signal or road load is described.

In FIG. 1, an illustrative example of a multi-channel road load input 2is shown. In the example shown, the multi-channel road load input 2 is adual channel input. The road load input 2 is time-varying, with timealong a first axis, e.g., an X-axis, and the channel magnitude along asecond axis, e.g., a Y-axis. However, skilled artisans should appreciatethat more than two channels, providing additional dimensions to the roadload input 2, and placed along different axes, may also be used withinthe scope of the disclosure.

The multi-channel road load input 2 is measured by sensors in thevehicle during service or computed as output from a vehicle dynamicscode. As nonlimiting examples, the sensors may be load sensors andtorque sensors. The sensors may be in wired or wireless communicationwith a data collection device (e.g., an external memory, standalonecomputer, networked computer, etc.) for later transmittal to the systemof the disclosure, or directly in wired or wireless communication withthe system of the present disclosure, as desired. One of ordinary skillin the art may select suitable sensors for providing the input 2, asdesired.

In the embodiment shown in FIG. 1, the multi-channel road load inputincludes a first channel 4 and a second channel 6. In the case of therubber bushing, the first channel 4 may be a radial load measured by theload sensor and the second channel may be a conical rotation measured bya torque sensor. It should be understood that any force, displacement,moment, rotation, etc., measured by other types of sensors, may also beembodied by the multi-channel road load input 2, as desired.

Referring now to FIG. 2, the multi-channel road load input 2 isprocessed into an alternate format for purposes of the analysisaccording to the method of the disclosure. In particular, the firstchannel 4 is disposed on a first axis, e.g., an X-axis, and the secondchannel 6 is disposed on a second axis, e.g., a Y-axis. Additionalchannels, if provided, are likewise placed on additional axes. It hasbeen surprisingly found that the arrangement of data from the pluralityof channels 4, 6 advantageously permits the fatigue analysis as shown inFIG. 3, and described further herein. It should also be appreciated thatthe arrangement of data from the plurality of channels 4, 6, as shownand described, provides certain advantages over rectangularly griddedconstructs, for which convergence difficulties at grid extremes may bean issue.

The method of the present disclosure, as shown in FIG. 4, includes astep 100 of providing the time-varying load data signal, for example, inthe form of the multi-channel road load input 2. Independent variablesof the time-varying load data signal are then identified. In particular,the independent variables may be the load signal and provided by thefirst channel 4 and the conical rotation signal provided by the secondchannel 6. Advantageously, the independent variables identified arevariables necessary for a subsequent damage calculation to determine apotential failure location in the elastomeric component. A skilledartisan may identify suitable independent variables for use in thedamage calculation, as desired.

Based at least in part on the identified independent variables, a step102 of conducting a finite element analysis (FEA) is performed in orderto obtain a base state 8 of an elastomeric component, for example, asshown in FIG. 3. The base state 8 of the elastomeric component ismodeled by FEA under simulated static or constant conditions, forexample, where the elastomeric component is assumed to be under a staticload. The FEA solution produces a database giving strain and stress atevery point or element in the model for the given base state. One ofordinary skill in the art understands how to conduct the FEA to obtainthe base state of the elastomeric component under static conditions.

In a next step 104, a plurality of case vectors 10, 12, 14, 16, 18 arethen selected. The case vectors 10, 12, 14, 16, 18 represent a space ofpossible loading states of the elastomeric component that may occurwithin the time-varying load data signal. Each of the case vectors 10,12, 14, 16, 18 begins at the base state 8, and radiate outwardlytherefrom.

Although five case vectors 10, 12, 14, 16, 18 are shown in FIG. 3 forpurposes of illustrating the method and system of the disclosure, itshould be understood that a fewer number or a greater number of the casevectors 10, 12, 14, 16, 18, in any desired range of orientations, mayalso be used within the scope of the disclosure. Furthermore, the casevectors 10, 12, 14, 16, 18 may be disposed along different dimensions,which in turn may be defined by additional channel inputs, as desired.

For at least a portion of the case vectors 10, 12, 14, 16, 18, a step106 is then performed in which a nonlinear FEA (as opposed to a “scaleand combine” or linear FEA) is conducted at discrete gridpoints 20, 22along the case vectors 10, 12, 14, 16, 18. The discrete gridpoints 20,22 may be selected at any location on the case vectors 10, 12, 14, 16,18 starting at the base state 8 and tracking the case vectors 10, 12,14, 16, 18. For example, the discrete gridpoints 20, 22 may bedistributed substantially evenly apart along a length of each of thecase vectors 10, 12, 14, 16, 18, randomly distributed, or selectivelydistributed for optimum analysis with respect to a particular desiredcurrent state 24.

In the example shown in FIG. 3, each of the discrete gridpoints 20, 22results in local solution variables (e.g., stresses, strains, nodaldisplacements, etc.) for every element from the FEAs conducted at eachof the discrete gridpoints 20, 22. Any suitable number of discretegridpoints 20, 22 may be selected along each of the case vectors 10, 12,14, 16, 18, as desired. It should be understood that the method andsystem may be further refined by adding additional case vectors 10, 12,14, 16, 18 together with modeling by FEA additional discrete gridpoints20, 22 along the case vectors 10, 12, 14, 16, 18.

Where the discrete gridpoints 20, 22 along the case vectors have beenmodeled by FEA, and local solution variables obtained from the FEA atthe discrete gridpoints 20, 22 in a step 108, a step 110 may then beemployed to obtain desired local solution variables (e.g., strain andstress) for the desired current state 24. The step 110 is performed byan interpolation engine 200 reading output from the FEA, for example, asshown in FIG. 6 and described further herein below. In particular, thedesired local solution variables permit an interpolation of strainhistory and stress history for the modeled elastomeric component in astep 110.

Referring now to FIG. 5, a suitable method for obtaining the desiredlocal solution variables for the desired current state 24 is shown insteps 110.1 to 110.5. For example, in step 110.1, the number of inputchannels 4, 6 is determined and assigned variable (N). In step 110.2,the desired current state 24 as a function of time is assigned a valueI_(i)(t), where (I) is the value of the input channel, (i) is thechannel, and time is the variable (t).

In step 110.3, the case vectors 10, 12 that neighbor the desired currentstate 24, I_(i)(t) are identified. As a nonlimiting example, theneighboring case vectors 10, 12 may be the nearest of the case vectors10, 12, 14, 16, 18 to the desired current state 24, I_(i)(t). However,it should be understood that certain case vectors 10, 12, 14, 16, 18 notnecessarily the nearest to the desired current state 24, I_(i)(t) mayalso be selected within the scope of the disclosure.

An interpolation cell 26, bounded by an upper and lower simplex definingthe interpolation cell 26, with edges that include the neighboring casevectors 10, 12, is then identified in step 110.4. It should beappreciated that the use of simplices permit a formation of theinterpolation cell 26 in any number of dimensions defined by any numberof channels. The identification of the interpolation cell 26 in step110.4 is conducted by identifying upper and lower discrete gridpoints20, 22 on the neighboring case vectors 10, 12, that bound the desiredcurrent state 24, I_(i)(t).

Where the interpolation cell 26 has been identified, it should beunderstood that the dependent variables at the desired current state 24,I_(i)(t) may be evaluated in a step 110.5, for example, using convenientinterpolation functions (e.g., piecewise constant interpolation, linearinterpolation, polynomial interpolation, spline interpolation, etc.).Once the interpolation cell 26 has been identified, other forms ofinterpolation including non-linear interpolation and combinations ofdifferent interpolation strategies may also be used, as desired. Thedesired current state 24, I_(i)(t) is thereby interpolated.

It should be appreciated that, in evaluating the dependent variables atthe desired current state 24, I_(i)(t), constraints such asincompressibility may also be enforced. Such constraints generally causethe results to comply with real-world limitations, and may be selectedby a skilled artisan, as desired.

In a step 112, a damage calculation based the interpolated strainhistory and/or the interpolated stress history to determine fatigue lifeand a potential failure location in the elastomeric component at thedesired current state 24 is then performed. In one example, the damagecalculation may be performed as described in U.S. Pat. No. 6,634,236 toMars, the entire disclosure of which is hereby incorporated herein byreference. Other means for performing the damage calculation using theinterpolated strain history and the interpolated stress history may alsobe employed within the scope of the present disclosure.

With renewed reference to FIG. 6, the interpolation engine 200 of thepresent disclosure is shown. The interpolation engine 200 is providedfor analysis of the time-varying load data signal obtained frommeasurement of the elastomeric component in service. The interpolationengine includes at least one processor 202 and at least one memory 204.The at least one memory 204 may be in the form of a tangible,non-transitory computer readable medium having a set ofcomputer-readable instructions 206 embodied thereon. The instructions206, when executed by the at least one processor 202, cause the at leastone processor 202 to perform the method of the present disclosure.

The memory 204 of the interpolation engine 200 includes at least oneinternal database. The internal database of the memory 204 has a uniquestructure. For each finite element, a base state (both independent anddependent variables) is stored in the database. For each case vector 10,12, 14, 16, 18 computed by the processor 202, the unit vector of thecase vector 10, 12, 14, 16, 18 is also stored in the database. For eachdiscrete case vector gridpoint 20, 22, independent variables (e.g., afirst parameter and a second parameter) and dependent variables (e.g.,strain tensor components and nodal displacements) are also stored in thedatabase.

In particular embodiments, the at least one processor 202 is configuredto obtain the desired local solution variables for each point or elementfrom the FEA of the elastomeric component based on the time-varying loaddata signal, i.e., the multi-channel road input 2 including the firstchannel 4 and the second channel 6, as described further hereinabove.

The at least one processor is also configured to provide at least one ofthe interpolated strain history and the interpolated stress history forfurther use in performing the damage calculation to determine thefatigue life and potential failure location in the elastomericcomponent. For example, as shown in FIG. 6, interpolated components ofnominal strain tensors (NE) as a function of time may be the output ofthe interpolation engine 200. It should be appreciated that theinterpolation of the system and method is particularly efficient, asopposed to running the finite element model through an entirety of theroad load input is prohibitively time-consuming.

The instructions 206 for execution by the at least one processor 202 mayalso be used to perform the other various steps of the method of thepresent disclosure, or permit the user to perform various steps, asdetailed further hereinabove.

The processor 202 may also be in communication with a human interface208, for example, at least one of a keyboard, a mouse, a video screen, atouch screen, and the like. The human interface 208 permits a user tointeract with the interpolation engine 200, for example, by providinginputs 210 for the selection of suitable case vectors 10, 12, 14, 16, 18and the creation of the plurality of FEA models 212 in accordance withthe disclosed method. The human interface 208 may also permit the userto upload the time-varying load data from measurements of theelastomeric component to the system. Other interactions between the userand the system as described hereinabove may be facilitated through useof the human interface 208.

In other embodiments, the time-varying load data may be uploaded to thesystem through a wired or wireless connection. For example, the systemmay be in communication with a computer network such as the Internet,through which the time varying load data is transmitted to the system.In another example, the time varying load data in the form of the firstchannel 4 and the second channel 6, saved from testing of theelastomeric component, may be uploaded to the system from a memorydevice such as a USB drive or the like. Other means for uploading thetime-varying load data to the system may also be used within the scopeof the present disclosure.

Advantageously, the method and system of the present disclosure providesa way to estimate local history of stress and strain, based on FEAmodeling of a series of load cases. It is adapted for the typical casewhere dynamic load perturbations are imposed on top of a static loadedstate. The method and system works with any number of independent inputchannels, and dependent local solution variables.

It is surprisingly found that the method and system of the disclosureavoids a need to model full time history in FEA. It accounts properlyfor material and kinematic nonlinear behavior, which is a feature notfound in other known failure analysis methods and systems. Moreover, thepresent method and system interpolates within a multidimensional space(e.g., one dimension for each input channel), which permitsinterpolation accuracy to be increased incrementally by addingadditional case vectors.

While certain representative embodiments and details have been shown forpurposes of illustrating the invention, it will be apparent to thoseskilled in the art that various changes may be made without departingfrom the scope of the disclosure, which is further described in thefollowing appended claims.

What is claimed is:
 1. A method for analyzing fatigue life of anelastomeric component, the method including steps operated by a systemincluding a user interface and an interpolation engine, theinterpolation engine in communication with the user interface, and theinterpolation engine including at least one processor and at least onememory having a tangible, non-transitory computer readable medium, thesteps comprising: receiving, by the interpolation engine, a base finiteelement analysis (FEA) with a base state of the elastomeric component,the base state of the elastomeric component obtained by analyzing theelastomeric component including at least one of stress and strain for atleast one element from the base FEA of the elastomeric component undersimulated static or constant conditions; generating, by theinterpolation engine, a predictive model including a multi-channelinput, the base state from the base FEA of the elastomeric componentunder simulated static or constant conditions, and generated casevectors to represent a space of loading states that occur within atime-varying load data signal of the multi-channel input, the casevectors originating at and extending outwardly from the base state inthe predictive model, the case vectors including a set of case vectorsneighboring a desired current state in the predictive model, thepredictive model further including an interpolation cell that bounds thedesired current state, the interpolation cell having edges that includethe set of the case vectors neighboring the desired current state,wherein the multi-channel input has independent variables; for the setof the case vectors neighboring the desired current state in thepredictive model, receiving by the interpolation engine from the userinterface a plurality of additional FEA of the elastomeric componentunder simulated dynamic conditions at a plurality of discrete gridpointsalong the set of the case vectors neighboring the desired current state,starting at the base state and tracking the set of the case vectorsneighboring the desired current state; obtaining, by the interpolationengine, local solution variables from the additional FEA of theelastomeric component at each of the plurality of discrete gridpoints,the local solution variables including at least one of stress and strainfor at least one element from the additional FEA for each of theplurality of discrete gridpoints, wherein the interpolation engine has adatabase included in the at least one memory that store the base FEAused to obtain the base state of the elastomeric component and each ofthe additional FEA used to obtain the local solution variables at eachof the discrete gridpoints; interpolating, by the interpolation engine,using the local solution variables associated with the interpolationcell and an interpolation function, at least one of a strain history anda stress history of the elastomeric component at the desired currentstate; and performing, by a set of instructions embodied on the at leastone memory of the interpolation engine and executed by the at least oneprocessor, a damage calculation based on one of the strain history andthe stress history interpolated for the elastomeric component, wherebythe fatigue life of the elastomeric component is predicted.
 2. Themethod of claim 1, further comprising a step of adding additional onesof the case vectors to increase interpolation accuracy of the predictivemodel.
 3. The method of claim 1, wherein the case vectors extend in twoor more dimensions from the base state.
 4. The method of claim 1,wherein the set of case vectors neighboring the desired current stateand interpolation weighting functions are evaluated and stored in thememory of the interpolation engine, prior to evaluation for each finiteelement, for purposes of enhanced execution speed.
 5. The method ofclaim 1, wherein the interpolation engine further generates a graphicalillustration of the predictive model at the user interface, thegraphical illustration including a grid with a first channel along afirst axis, and a second channel along a second axis, the base stateidentified on the grid and having the generated case vectors radiatingoutwardly therefrom, the discrete gridpoints shown on the set of thecase vectors neighboring the desired current state, and the desiredcurrent state shown disposed within the interpolation cell having theedges that include the set of the case vectors neighboring the desiredcurrent state.
 6. The method of claim 1, further comprising a step ofmaking another elastomeric component with a parameter modified inresponse to the fatigue life prediction.
 7. A system for analyzingfatigue life of an elastomeric component, comprising: a user interface,and an interpolation engine in communication with the user interface,and the interpolation engine including at least one processor and atleast one memory having a tangible, non-transitory computer readablemedium, the interpolation engine performing the steps comprising:receiving, by the interpolation engine, a base finite element analysis(FEA) with a base state of the elastomeric component, the base state ofthe elastomeric component obtained by analyzing the elastomericcomponent including at least one of stress and strain for at least oneelement from the base FEA of the elastomeric component under simulatedstatic or constant conditions; generating, by the interpolation engine,a predictive model including the multi-channel input, the base statefrom the base FEA of the elastomeric component under simulated static orconstant conditions, and generated case vectors to represent a space ofloading states that occur within the multi-channel input, the casevectors originating at and extending outwardly from the base state inthe predictive model, the case vectors including a set of case vectorsneighboring a desired current state in the predictive model, thepredictive model further including an interpolation cell that bounds thedesired current state, the interpolation cell having edges that includethe set of the case vectors neighboring the desired current state; andfor the set of the case vectors neighboring the desired current state inthe predictive model, receiving by the interpolation engine from theuser interface a plurality of additional FEA of the elastomericcomponent under simulated dynamic conditions at a plurality of discretegridpoints along the set of the case vectors neighboring the desiredcurrent state, starting at the base state and tracking the set of thecase vectors neighboring the desired current state; obtaining, by theinterpolation engine, local solution variables from the additional FEAof the elastomeric component at each of the plurality of discretegridpoints, the local solution variables including at least one ofstress and strain for at least one element from the additional FEA foreach of the plurality of discrete gridpoints, wherein the interpolationengine has a database included in the at least one memory that store thebase FEA used to obtain the base state of the elastomeric component andeach of the additional FEA used to obtain the local solution variablesat each of the discrete gridpoints; interpolating, by the interpolationengine, using the local solution variables associated with theinterpolation cell and an interpolation function, at least one of astrain history and a stress history of the elastomeric component at thedesired current state; and performing, by a set of instructions embodiedon the at least one memory of the interpolation engine and executed bythe at least one processor, a damage calculation based on one of thestrain history and the stress history interpolated for the elastomericcomponent, whereby the fatigue life of the elastomeric component ispredicted.
 8. The system of claim 7, wherein the interpolation enginefurther generates a graphical illustration of the predictive model atthe user interface, the graphical illustration including a grid with afirst channel along a first axis, and a second channel along a secondaxis, the base state identified on the grid and having the generatedcase vectors radiating outwardly therefrom, the discrete gridpointsshown on the set of the case vectors neighboring the desired currentstate, and the desired current state shown disposed within theinterpolation cell having the edges that include the set of the casevectors neighboring the desired current state.
 9. The system of claim 7,further comprising the step of making another elastomeric component witha parameter modified in response to the fatigue life prediction.
 10. Thesystem of claim 7, wherein the case vectors extend in two or moredimensions from the base state.